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Q , Abstract 

Q ■ We consider a D7-brane probe of AdSs x S in the presence of pure gauge 5-field. 

. In the dual gauge theory, the -B-field couples to the fundamental matter introduced 

' by the D7-brane and acts as an external magnetic field. The S-field supports a 

6-form Ramond-Ramond potential on the D7-branes world volume that breaks the 
supersymmetry and enables the dual gauge theory to develop a non-zero fermionic 
condensate. We explore the dependence of the fermionic condensate on the bare quark 
' mass mq and show that at zero bare quark mass a chiral symmetry is spontaneously 

broken. A study of the meson spectrum reveals a coupling between the vector and 
^ i scalar modes, and in the limit of weak magnetic field we observe Zeeman splitting of 

^ \ the states. We also observe the characteristic ^Jfn^ dependence of the ground state 

corresponding to the Goldstone boson of spontaneously broken chiral symmetry. 
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1 Introduction 



In recent years, progress has been made in the study of gauge theory with matter in 
the fundamental representation in the context of gauge/string duahties generahzing the 
AdS/CFT correspondence. One way to achieve this is by introducing D7-branes in the 
probe hmit[2j that amounts to the condition Nf -C Nc- The fundamental strings stretched 
between the stack of Nc D3-branes and the Nf flavor D7-branes give rise to M'=2 hypermul- 
tiplet. The separation of the D3- and D7-branes in the transverse directions corresponds 
to the mass of the hypermultiplet, the classical shape of the D7-brane encodes the value 
of the fermionic condensate, and its fluctuations describe the light meson spectrum of the 
theory [3]. This technique for introducing fundamental matter has been widely employed 
in different backgrounds. Of particular interest is the study of non supersymmetric back- 
grounds and phenomena such as spontaneous chiral symmetry breaking. These phenomena 
were first studied in this context in ref. [3], using analytical and numerical techniques. In 
several works this approach was further developed, and has proven itself a powerful tool for 
the exploration of gauge theories, in particular, for the description of their thermodynamic 
properties or for the building of phenomenological models relevant to QCD [5] [39]. 

In this paper we will be interested in introducing fundamental matter into the gauge 
theory in the presence of an external electromagnetic field that couples to the fundamental 
fermions. The supersymmetry will be explicitly broken by the external field, and we will 
observe spontaneous symmetry breaking, and non-trivial mixing in the spectrum of mesons. 



2 Fundamental Matter in an External Magnetic Field 
2.1 Basic Configuration 

Let us consider the AdSs x 5*^ geometry describing the near-horizon physics of a collection 
of Nc extremal D3-branes. 

ds"^ = ^{-dxl + dxj + dxl + dxj) + R"^^ + R'^dnj , (1) 
gsC(4) = ^^2;° A dx^ A dx^ A dx^ , 



e = gs , 
i?^ = 4:7ig,Nca 



Where cifig is the unit metric on a round 5^. In order to introduce fundamental matter we 
first rewrite the metric in the following form, with cifig the metric on a unit S*^: 

ds^ = f-—^ — [-dxl + dxl + dxl + dxl] + — — [dp^ + p^d^l + dL'^ + , 

R p -\- Li 

dnj = d^^ + cos'^tpdp'^ + sm'^^dj^, (2) 
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where p, ip^ (3, 7 and L, are polar coordinates in the transverse and respectively. 
Note that: v? = + L^. We use xo,xi,X2, X3, p,ip, P, 'J to parametrise the world volume 
of the D7-brane and consider the following ansatz [3] for its embedding: 

= const, L = L{p), 

leading to the following form of the induced metric on its worldvolume: 

ds = Pl±^PjL[-dxl + dxl + dxl + dx^] + -^-^^[(1 + L'{pf)dp' + pHnl] . (3) 

Now let us consider the general DBI action: 

Sdbi = -P7 j rf«ee~*[-det(G,fe + Bab + 2T:a'Fab)]"^ . (4) 

Ms, 

Here p^ = [(27r)''a'^]~^ is the D7-brane tension, Gab and Bab are the induced metric 
and -B-field on the D7-brane's world volume, while Fab is its world-volume gauge field. A 
simple way to introduce magnetic field would be to consider a pure gauge 5-field along 
parts of the D3-branes' world volume, e.g.: 

5^2) = Hdx2 A dxs . (5) 

Since Bab can be mixed with the gauge field strength Fab, this is equivalent to a magnetic 
field on the world-volume. Recently a similar approach was used to study drag force in 
SYM plasma [35J. Note that since the 5-field is pure gauge, dB = 0, the corresponding 
background is still a solution to the supergravity equations of motion. On the other hand, 
the gauge field Fab comes at next order in the a' expansion compared to the metric and 
the i?-field components. Therefore to study the classical embedding of the D-brane one 
can study only the {Gab + Bab) part of the DBI-action. However, because of the presence 
of the -B-field, there will be terms ot first order in a' in the full action linear in the gauge 
field Fab- Hence integrating out Fab will result in a constraint for the classical embedding 
of the D7-brane. 

Since for our configuration, we have that: 

^(2) ^ q{2) ^ Q ^ q{2) ^ fj^^^ ^ Q ^ 

and at first order in a' the only contribution to the Wess-Zummino is 

27ra'/i7 j F^2) A C(6) • (6) 
By using the following expansion in the DBI action: 

[-detiEab + 2WF„f,)]i/2 = na'VEE'"'Fab + 0{F^)- E = -deiEab] , (7) 
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where we have introduced Eat = Gab + Bab as a notation for the generahzed induced metric, 
we obtain the following action to first order in a': 



= 7ia'^ j d^^VEE^'^'^F^ab] + 27raV7 J F^^) A C(6) . (8) 

Ms 

The resulting equation of motion does not contain Aa and sets the following constraint for 
the C(6) potential induced by the gauge S-field. 

|^a6Ai...A6 5^C'^^ ,^ = -a,(v^E[H); a,b,flu . . . fie e Ms; . (9) 

Note that C(6) has a dynamical term proportional to I/kq in the supergravity action, and 
that the D7-brane action is proportional to /iy = 211/ Kq. Therefore they are at the same 
order in a' [H]. We must solve for C(6) using the action: 



(10) 



The solution obtained from equation (fTOj) has to satisfy the constraint given in equation ([9]). 
Our next goal will be to find a consistent ansatz for C(6). To do this let us consider the 
classical contribution to the DBl action: 

Sns = f d'^^/E . (11) 



From equation ( fTTl) one can solve for the classical embedding of the D7-brane, which 
amounts to second order differential equation for L{p) with some appropriate solution 
Lq{p). After substituting Lo(p) in (fTTl) we can extract the form of the C(6) potential 
induced by the S- field. However one still has to satisfy the constraint ([9]). It can be 
verified that with the choice ^ for the S-field and the ansatz of equation ^ for the 
induced metric, the right-hand side of equation ^ is zero. Then equation ^ and the 
effective action ffTOj) boil down to finding a consistent ansatz for C(6) satisfying: 



d,{V^dCf''^^^) = -^H5{L - L,{p)) , (12) 

TT 

^^(v^^iCr^^"'^) =-L',{p)^H5{L-L,{p)) , (13) 

TT 

gabMi...A6 5^C'^^ ,^ =0; a,6,/ii,.../i6 G^s; • (14) 



One can verify that the choice: 



C'(6)oipV>a/3 = \fiP^ ^' V^)' ^C'(6)^o,^^„^ = OlJ , (15) 
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is a consistent ansatz and the solution for the C(6) field strength can be found to be: 



^C'(6)iow = "T^^ L(p2 + L2)2 ®^^ " sin COS 7/- . (16) 

It is this potential which breaks the supersymmety. It is important to note that there 
is no contradiction between the fact that the 5-field that we have chosen does not break 
the supersymmetry of the AdSs x 5*^ supergravity background, on the one hand, and the 
fact that the physics of the D7-brane probing that background does have supersymmetry 
broken by the i?-field, on the other. This is because the physics of the probe does not 
back-react on the geometry. 

In what follows, we will study the physics of the D7-branes and the resulting dual gauge 
theory physics. Among the solutions for the D7-brane embedding, there will be a class 
with non-trivial profile having zero asymptotic separation between the D3- and D7-branes. 
This corresponds to a non-zero fermionic condensate at zero bare quark mass. Therefore 
the non-zero background magnetic field will spontaneously break the chiral symmetry. 
Geometrically this corresponds to breaking of the S0{2) rotational symmetry in the (L, 0)- 
plane [5]. 

2.2 Properties of the Solution 

We now proceed with the exploration of the properties of the classical D7-brane embedding. 
If we consider the action flTT]) at leading order in a', we get the following effective lagrangian: 



C = -^p^sm^cos^./l^Jl+ f . (17) 

9s V (p^ + L^y 

The equation of motion for the profile Lq{p) of the D7-brane is given by: 



As expected for large (Lq + p^) ^ oo or if — > 0, we get the equation for the pure AdSs x 5*^ 
background [2]: 

a„ I ^ ^0 1=0. 



^0 

Therefore the solutions to equation f|T8l) have the following behavior at infinity: 

Lo(p) =m + 4 + ---- (19) 
P 
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where the parameters m (the asymptotic separation of the D7- and D3- branes) and c (the 
degree of bending of the D7-brane) are related to the bare quark mass ruq = m/2TTa' and 
the fermionic condensate {^ip) oc — c respectively [5j. As we shall see below, the presence 
of the external magnetic field and its effect on the dual SYM provide a non vanishing value 
for the fermionic condensate, furthermore the theory exhibits chiral symmetry breaking. 

Now notice that H enters in f ll7l) only through the combination H^R^. The other natu- 
ral scale is the asymptotic separation m. It turns out that different physical configurations 
can be studied in terms of the ratio = m^/ [HR"^): Once the m dependence of our solu- 
tions are known, the m and H dependence follows. Indeed let us introduce dimensionless 
variables via: 

p = RVHp , Lo = RVhL , L[,(p) = L'{p) . (20) 
The equation of motion ( |T8l) then takes the form: 



The solutions for L(p) can be expanded again to: 

L{p)=m + ^ + ..., (22) 

and using the transformation fl20|) we can get: 

c = ~cR^H'/^ . (23) 

It is instructive to study first the properties of fl2T|) for m ^ 1, which corresponds to 
weak magnetic field H <^ m^/i?^, or equivalently large quark mass m ^ R\fH. 

2.2.1 Weak Magnetic Field 

In order to analyze the case of weak magnetic field let us expand L{p) = rh + r]{p) and 
linearize equation fl2T|) while leaving only the leading terms in (p^ + m^)~^. The result is: 



which has the general solution: 

Ci rh 
p2 4:p^{fh? + p 



= l^-I157SJ^ + '^- (25) 



Prom the definition of r]{p) and equation fl22l) we can see that Ci = c and since r]\ 
we have C2 = 0. Now if we consider rh large enough, equation ( !25l) should be valid for all 



p—tOO 
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p. It turns out that if we require that our solution be finite as p — we can determine the 
large fh behavior of c. Indeed the second term in (!25l) has the expansion: 

'fh 1 1 1 ^. o, , , 

+ + (26) 



Ap^im? + p^) Am p^ Am^ 
Therefore we deduce that: 



and finally, we get for the profile of the D7-brane for m ^ 1: 

Hp)=m+J^^- ,~2t~2^~2^ • (28) 

Amp'' Ap''[m'^ + p^) 

If we go back to dimensionful parameters we can see, using equations (1251) and flTTI) that 
for weak magnetic field H the theory has developed a fermionic condensate: 

(^^) oc -c = -^H' . (29) 
Am 

However this formula is valid only for sufficiently large m and we cannot make any 
prediction for the value of the fermionic condensate at zero quark mass. To go further, the 
involved form of equation fl2T|) suggests the use of numerical techniques. 



2.2.2 Numerical Results 

In this subsection we solve numerically equation fl2T|) for the embedding of the D7-brane, 
using Mathematica. It is convenient to use initial conditions in the IR as has been recently 
discussed in the literature [B], [T|.We use the boundary condition L'(p)|p=o = 0. We used 
shooting techniques to generate the embedding of the D7 for a wide range of fh. Having 
done so we expanded numerically the solutions for L(p) as in equation (l22l) and generated 
the points in the (m, — c) plane corresponding to the solutions. The resulting plot is 
presented in figure [H 

As one can see there is a non zero fermionic condensate for zero bare quark mass, the 
corresponding value of the condensate is = 0.226. It is also evident that the analytical 
expression for the condensate fl27|) that we got in the previous section is valid for large m, 
as expected. Now using equation fl23|) we can deduce the dependence of Ccr on H: 

Ccr = CcrR^H^''^ = 0.22QR^H^/^ . (30) 

It is interesting to check the consistency of our numerical analysis by solving equa- 
tion f|T8|) numerically and extracting the value of Ccr for wide range of R^H, the resulting 
plot fitted with equation fl30|) is presented in figure [2l 

Another interesting feature of our phase diagram is the spiral behavior near the origin 
of the (m, — c)-plane which can be seen in figure [31 A similar feature has been observed in 
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-c(m) 




m 



12 3 4 




Figure 1: The black line corresponds to (j27p . one can observe that the analytic result is valid for 
large m. It is also evident that for m = {ipip) ^ 0. The corresponding value of the condensate 
is Ccr = 0.226 



ref. [7j, where the authors have argued that only the lowest branch of the spiral correspond- 
ing to positive values of m is the stable one (corresponding to the lowest energy state). 
The spiral behavior near the origin signals instability of the embedding corresponding to 
Lq = 0. If we trace the curve of the diagram in figure [3] starting from large m, as we go 
to smaller values of m we will reach zero bare quark mass for some large negative value of 
the fermionic condensate Ccr- Now if we continue tracing along the diagram one can verify 
numerically that all other points correspond to embeddings of the D7-brane which inter- 
sect the D3-brane at least once. (Note also that in ref. [4j, such behavior was considered 
inconsistent with the interpretation of the embedding as a re- normalization group flow.) 
After further study one finds that the part of the diagram corresponding to negative values 
of m represents solutions for the D7-brane embedding which intersect the D3-branes odd 
number of times, while the positive part of the spiral represents solutions which intersect 
the D3-branes even number of times. The lowest positive branch corresponds to solutions 
which don't intersect the D3-branes and is the stable one, while the upper branches have 
correspondingly 2,4, etc., intersection points and are ruled out. 

3 Meson Spectrum 
3.1 General Properties 

We study the scalar meson spectrum. To do so we will consider quadratic fluctuations [3] 
of the embedding of the D7-brane in the transverse (L, </))-plane. It can be shown that 
because of the diagonal form of the metric the fluctuation modes along the (p coordinate 
decouple from the one along L. However, because of the non-commutativity introduced 
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Figure 2: A plot of the magnitude of the fermionic condensate at zero bare quark mass Ccr as 
function of R^H, the black curve represents equation (j30p . 
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Figure 3: A magnification of figure [T] to show the spiral behavior near the origin of the (— c, m)- 
plane. 



by the -B-field we may expect the scalar fluctuations to couple to the vector fluctuations. 
This has been observed in ref. [8], where the authors considered the geometric dual to non- 
commutative super Yang Mills . In our case the mixing will be even stronger, because of the 
non-trivial profile for the D7-brane embedding, resulting from the broken supersymmetry. 

Let's proceed with obtaining the action for the fluctuations. To obtain the contribution 
from the DBI part of the action we consider the expansion: 

L = Lo{p) + 27ra'x, </> = + 27ra' , (31) 

where Lo(p) is the classical embedding of the D7-brane solution to equation f[T8l) . To second 
order in a' we have the following expression: 

E,, = El + 2ixa'El, + {2ixa'fEl , (32) 
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where E , E , E are given by: 



=Gab{p,Lo{p),^) + B^, 

R^Lq 



p2 _ 



2 , t2 i^aXK + 9bXK) + dloGabX + Fab 
P + ^0 



(33) 



{daxK + dbX5'a)x+ldlGabx' 



Here Gab and Sab are the induced metric and B field on the D7-brane's world volume. Now 
we can substitute equation (!33l) into equation (fTTj) and expand to second order in a'. It is 
convenient [8J to introduce the following matrices: 



\Eab\ 



-' = s + j, 

where S is diagonal and J is antisymmetric: 

Gil 



jab 

Gn 

G rvrv 



diag{-Gn , 



Gn 



Gl, + H^' Gl, + H^' 



5 '^pp 5 '-'i/ji/;' '-^aof' ^0l3i 



f3(3S ' 



^11 



2 - '53'^2) , 



pp 



^2(1 + ^0 



G 



tpip 



r2' 
-^0 



cos^ ipG^^ ; Gfsp = sin^ i>G^^ . 



(34) 

(35) 
(36) 

(37) 



Now it is straightforward to get the effective action. At first order in a' the action for 
the scalar fluctuations is the flrst variation of the classical action fill I) and is satisfled by 
the classical equations of motion. The equation of motion for the gauge fleld at flrst order 
was considered in Section 2 for the computation of the G(6) potential induced by the B- 
fleld. Therefore we focus on the second order contribution from the DBI action. 

After integrating by parts and taking advantage of the Bianchi identities for the gauge 
fleld, we end up with the following terms. For x'- 



and for F: 



S 



ab 



p2 + Lo'l + L: 



j2^aXdbX + 



L' 



^X^(3. 



1 



and the mixed terms: 



Cf oc -V^S'^'^'S^'^'FabFa^y 



^Fx a: — 7, — /X^23 



(39) 



(40) 
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and for $: ^ 

oc l^^^I^S'^'da'^db'^ , (41) 



where the function / in (HOl) is given by: 

l + Lo 



f{p) = d, 9{p)—^J'' ] + J''d:^o9ip) + 2g{p)j''S''dL,,G,, , (42) 



with g{p) =-^^ - = pVl + + 



sin cos -i/^ y (P^ + -^o)^ 

As can be seen from equation ( llOi) the ^2,^3 components of the gauge field couple 
to the scalar field x via the function /. Note that since for p — 00 and L — >■ 00, we 
see that J^^ — 0, the mixing of the scalar and vector field decouples asymptoticly. In 
order to proceed with the analysis we need to take into account the contribution from the 
Wess-Zumino part of the action. The relevant terms to second order in a' arc [8j: 

Swz = -^^^/i? J i^(2) A F(2) A ^4) + (27r«')/i7 J i^(2) A B^^) A P[C^,)] , (43) 

where C(4) is the background R-R potential given in equation ([1]) and C(4) is the pull back 
of its magnetic dual. One can show that: 

= — 7-^ ^-^^ sin ip cos i/jdi/j A da A d/S A d(f) . (44) 

gs [p^ + L^Y 

Writing = 27rQ;'$ we write for the pull back P[C(4)]: 

P[q4)] = ^K{p)daMi) A da A dp A dx", (45) 



K{P) = -R'Ll^^^ (46) 



where we have defined: 

34 r2 2p + Lq^ 

Now note that the S-field has components only along and x^, therefore dx'^ in equation 
( l45ll can be only dp,dx^ or (ix^. This will determine the components of the gauge field 
which can mix with $, However after integrating by parts and using the Bianchi identities 
one can get the following simple expression for the mixing term: 

- (2W)^^ / d'^'^Hd,K<^Fo, , (47) 

resulting in the following contribution to the complete lagrangian: 

Lp^ oc ^^/79,K<l>Foi . (48) 
10 



Note that this means that only the Aq and Ai components of the gauge field couple to the 
scalar field $. Next the contribution from the first term in (143|) is given by: 



(2vra')^^ / d'^ ^P F^,F^,e^'-^ , (49) 

gs J 

where the indices take values along the p,ip,a,f3 directions of the world volume. This 
will contribute to the equation of motion for Ap, A^, Aa and Af^, which do not couple to 
the scalar fiuctuations. In this paper we will be interested in analyzing the spectrum of 
the scalar modes, therefore we will not be interested in the components of the gauge field 
transverse to the D3-branes world volume. However although there are no sources for 
these components from the scalar fiuctuations, they still couple to the components along 
the D3-branes as a result setting them to zero will impose constraints on the Aq . . . A^. 
Indeed from the equation of motion for the gauge field along the transverse direction one 
gets: 

3 

J2S''''dtdaAa = 0, b = p,^,a,f3, (50) 

a=0 

(Here, no summation on repeated indices is intended.) However the non-zero i?-field ex- 
plicitly breaks the Lorentz symmetry along the D3-branes' world volume. In particular we 
have: 

which suggests that we should impose: 

- ^0^0 + diA^ = , ^2^2 + = . (52) 

We will see that these constraints are consistent with the equations of motion for Aq . . . A^. 
Indeed with this constraint the equations of motion for x, ^ and yl^, /i = . . . 3 are, for x'- 

''''\(.^)^^^r^-x. (53) 



1 + f ^ f dg L^, \ d'g\ , 1 + L', 



g V Kdlol + L'^'J dLf^J g 

and for $: 

and finally for Aa- 

1 ^ ^ -^D^o + ^a,* ^ , (55) 
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--dsx = 
9 



= . 

9 



We have defined: 



92 + 92 



^ = -dt + di + -^-^. (56) 

As one can see the spectrum sphts into two independent components, namely the vector 
modes Aq^Ai couple to the scalar fiuctuations along $, while the vector modes v42,v43 
couple to the scalar modes along x- However it is possible to further simplify the equations 
of motion for the gauge field. Focusing on the equations of motion for Aq and Ai in 
equnation ( l55l) . it is possible to rewrite them as: 



^ V 1 + ^0 / P (P + ^0)' 

Note that the first constraint in (l52l) trivially satisfies the second equation in (!57j) . In this 
way we are left with the first equation in (1571) . Similarly one can show that using the 
second constraint in fl52l) the equations of motion in (|55|) for A2 and A3 boil down to a 
single equation for F23: 

9dpF23 \ 




:i+^o^)(i + (^)y (P^+^g)^+M^ 

+ 2nt"^^'4^ +kdl + dl)x = 0. (58) 
Now let us proceed with a study of the fiuctuations along $. 



3.2 Fluctuations Along $ 

To proceed, we have to take into account the Fqi component of the gauge field strength 
and solve the coupled equations of motion. Since the classical solution for the embedding 
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of the D-brane is known only numerically we have to rely again on numerics to study the 
meson spectrum. However if we look at equation (fT8|) we can see that the terms responsible 
for the non-trivial parts of the equation of motion are of order H^. On the other hand, 
the mixing of the scalar and vector modes due to the term (HSj) appear at first order in H. 
Therefore it is possible to extract some non-trivial properties of the meson spectrum even 
at linear order in H and as it turns out, we can observe a Zeeman-like effect: A splitting 
of states that is proportional to the magnitude of the magnetic field. To describe this, let 
us study the approximation of weak magnetic field. 

3.2.1 Weak magnetic field 

To first order in H the classical solution for the D7-brane profile is given by: 

Lo = m + 0{H^), (59) 

where m is the asymptotic separation of the D3 and D7-branes and corresponds to the 
bare quark mass. In this approximation the expressions for g{p) and dpK{p), become: 

g{p) = , dpK{p) 



and the equations of motion for $ and Fqi, equations (1541) and simplify to: 

- pW9p$ + JlL^^ □$-4g 2^3 -^01=0, 60 

and -^d, {p%F,,) + + , ,5 2,2 °-^oi - fJ^\,, V'^ = , 

p-^ ^ p^ [p + m'^Y [p^ + m^y 

where □ = -9^ + dl + dl + (9|, = -dl + dl . 

This system has become similar to the system studied in ref. [8] and in order to decouple 
it we can define the fields: 

0± = Foi ± mV^ , (61) 



where V = \/—dQ + d\. The resulting equations of motion are: 

-d,{p%<j)±) + + , . ,., a0± T H V<P± = . (62) 

p-^ l^pA _|_ yji'^Y l^p^ _|_ m^'j'i 

Note that is the Casimir operator in the {xo,Xi) plane only, while □ is the Casimir 
operator along the D3-branes' world volume. If we consider a plane wave e*^ '^ then we can 
define: 

and we have the relation: 

= Mo\ -kl-kj . (64) 
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The corresponding spectrum of is continuous in k2.,k^. However, if we restrict ourselves 
to motion in the (xq, a;i)-plane the spectrum is discrete. Indeed let us consider the ansatz: 

0± = r7±(p)e-'"°'=°+*^i"i . (65) 

Then we can write: 

j^d,{p%r,^) + ^^,f^,y Mir^± T H^0^^M^r^^ = , (66) 
M± = Moi± . 

Let us analyze equation (!66|) . It is convenient to introduce: 

2/ = -4; M± = —M^- P±{y) = {l-yT^v±. (67) 



R 



2a± = 1 + -v/l + Ml; e = H—. 



With this change of variables equation (!66l) is equivalent to: 

y{l - y)P'i + 2(1 - (1 - a±)y)P' - a±{a±.,)P± ± e-^^^^P± = . (68) 

Next we can expand: 

P± = Po ± ePi + 0{e^) ; a± = ± eai + 0{e^) ; (69) 
M± = Mo±ea^^^^^^ + 0{e') ; Mq = 2^ao{ao + 1) . 

IVIo 

leading to the following equations for Pq and Pi: 

y{l-y)P;; + 2{l-{l-ao)y)P;,-ao{ao-l)Po = 0, (70) 
yil-y)P[' + 2{l-{l-ao)y)P[-ao{ao-l)Pi = (ai(2ao - 1) 



The first equation in fl7Up is the hypergeometric equation and corresponds to the fluctua- 
tions in pure AdS^ x S'^ . It has the regular solution [3]: 

Poiy) = F{-ao,l-ao,2,y) . (71) 

Furthermore regularity of the solution for r]{p) at infinity requires [3J that ao be discrete, 
and hence the spectrum of Mq: 

1 — ao = — n, n = 0, 1, . . . (72) 
Mo = 2v/(n + l)(n + 2) . 
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The second equation in (!70|) is an inhomogeneous hypergeometric equation. However for 
the ground state, namely n = 0, Pq = F{—1, 0, 2,y) = 1 and one can easily get the solution: 



Pi{y) = ^ ln(l -y) + {Qa, - Mo)(ln(-2/) + ^)- 



(73) 



On the other hand, using the definition of P±{y) in fl67|l to first order in e we can write: 

= (l T ln(l - y)) (1 ± ePM) , (74) 

for the ground state ao = 1 and we end up with the following expression for rj±: 



±e 



1-y 4(1 -y)2 1-y 



y 



6 



(75) 



Now if we require that our solution is regular at y = and goes as 1/p^ oc l/y at infinity, 
the last term in (1751) must vanish. Therefore we have: 



Mo 
6 



(76) 



After substituting in (169!) and (1671) we end up with the following correction to the ground 



M 



10 
8 
6 
4 



n=3 



n=2 



n=0 



0.05 



0.1 



0.15 



0.2 



H/m 



Figure 4: Plot of M = MR'^/m vs. iJ/m for the first three states. The dashed black lines 
correspond to the spectrum given by equation ([72]) 



sate: 



M± = Mo ± 



H 



m 



(77) 
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We observe how the introduction of an external magnetic field breaks the degeneracy of 
the spectrum given by equation (!72|) and results in Zeeman splitting of the energy states, 
proportional to the magnitude of H. Although equation (177j) was derived using the ground 
state it is natural to expect that the same effect takes place for higher excited states. To 
demonstrate this it is more convenient to employ numerical techniques for solving equation 
and use the methods described in ref. pE] to extract the spectrum. The resulting plot 
is presented in figure HI As expected we observe Zeeman splitting of the higher excited 
states. It is interesting that equation (1771) describes well not only the ground state, but 
also the first several excited states. 

It turns ou that one can easily generalize equation (1771) to the case of non-zero mo- 
mentum in the {x2-, X3)-plane. Indeed if we start from equation ( l62l) and proceed with the 
ansatz: 

0± = r/±(p)e-^^-'= , (78) 



we end up with: 



1 AR^Tn 

y^d,{p%v^) + yr^^.Mlf,^ T Hyr^^,M,,^V± = , (79) 
Moi+ = Jul + H,: A;23 = xkl + B . 



After going through the steps described in equations (167 1) -( 1751) . equation ( 1761) gets modified 
to: 



Mo L, kj. 



Note that validity of the perturbative analysis suggests that ai is of the order of ao and 
therefore we can trust the above expression as long as /C23 is of the order of Mq. Now it is 
straightforward to obtain the correction to the spectrum: 



M4 



We see that the addition of momentum along the {x2 — X3)-plane enhances the splitting of 
the states. Furthermore the spectrum depends continuously on ^23. 




3.2.2 Strong Magnetic Field 

For strong magnetic field we have to take into account terms of order H^, which means 
that we no longer have an expression for Lq{p) in a closed form and we have to rely on 
numerical calculations only. Furthermore there is no obvious way to decouple equations 
( I54p and (IFT]) . However, it is still possible to extract information about the spectrum of 
the scalar modes if we restrict ourselves to fiuctuations along the {x'^,x^) plane. In this 
way there is no source term in equation (l55l) . and we can consistently set Fqi equal to 
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zero. We consider time independent fluctuations satisfying tlie ansatz e""'^^''^^, (wliere 
is tlie radial coordinate in the {x2 — X3)-plane). The damping factor in the exponent can 
be thought of as the mass of the scalar meson in 2 Euclidean dimensions. Indeed let us 
consider the ansatz: 

$ = /i(p)e-^^2^'-^'=3^Y;(S^) , (82) 

where Yi{S^) are the spherical harmonics on the 5*^ sphere satisfying: A^gY; = + 2)Yi. 
With this set-up the equation of motion for $, equation fl5^ . reduces to equation for h{p): 



Ittit ) — ^ + {p- + Lif + R^m^^p^ - ° ' 

where we have defined: 

^^23 = -kl - kl . (84) 

Before we proceed with the numerical analysis of equation fl55]) let us introduce dimension- 
less variables by performing the transformation (l20l) and defining: 

= ^^^23 . (85) 

The resulting equation is: 



H 



p3^1 + L''(l + (p2 + L2)2)l/2 V ^ (p2 + I'')V V/TTX^ 

- ^/(/ + 2)Mp)+ -^'^g ^ Mp)=0- (86) 

(p2 + L2)2 + 1 

In order to study the spectrum we look for normalizable solutions which have asymptotic 
behavior /i(p) oc 1/p^ for large p and satisfy the following boundary conditions at p = 0: 

/i'(0) = 0; /i(0) = 1 . (87) 

Let us consider first the lowest level of the spectrum. The spectrum that we get as a 
function of the bare quark mass is plotted in figure [51 

For large m the spectrum asymptotes (the dashed line in figure [5]) to the one for pure 
AdS'^ X space obtained by ref. [3] 

Mo = ^V(r2 + / + l)(n + / + 2) , (88) 

with the substitution n = 0, / = 0, to obtain our case. Therefore we are describing the 
lowest possible state of the meson spectrum. In figure [6] we have zoomed in the area near 
the origin of the (m, m23)-plane, one can see that for small values of m = 2Txa'mq/ R\/H 
we observe oc ^Jrf^ dependence of the ground state on the bare quark mass m^, which is to 
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Figure 5: Spectrum of 17123 vs. rh. The dashed hne represents the lowest level of the meson 
spectrum for pure AdSs x space 




0.0025 0.005 0.0075 0.01 0.0125 0.015 0.0175 



Figure 6: Enlargement of part of the spectrum of vs m from figure O The black solid curve 
shows the oc \/m fit. 



be expected since the chiral symmetry associated with the spinor representation of SO (2) 
is spontaneously broken [ID]. 

It is interesting to look for modes corresponding to higher excited states (non-zero n). 
In figure [7] we have presented a plot of some of these. Again, the dashed line correspond to 
the pure AdS^ x 5*^ spectrum given by flHHl) for / = 0. For small values of m one can observe 
the qualitative difference of the behavior of the spectrum corresponding to the n = state 
from that of the higher excited states. Indeed as m — > the n = states follow the y/^ 
behavior plotted in figure El while the excited states tend to some finite values at zero bare 
quark mass. The n = states merge into the Goldstone boson of the spontaneously broken 
chiral symmetry. 
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Figure 7: Spectrum of 77123 vs rh for n = ... 4. The dashed lines represent the spectrum for 
AdS^ X space. 
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